Neutron Stars OS Qj Probe for Dense Matter 



V.A. Dexheimer a S. Schramm a ' b 

a FIAS, Johann Wolfgang Goethe University, Frankfurt am Main, Germany 
h CSC, Johann Wolfgang Goethe University, Frankfurt am Main, Germany 



Abstract 

We study different stages of the neutron star cooling by computing neutron star properties at 
various temperatures and entropies using an effective chiral model including hadronic and quark 
degrees of freedom. Macroscopic properties of the star such as mass and radius are calculated 
and compared with observations. It can be seen that the effects of chiral restoration and de- 
confinement to quark matter in the core of the neutron star at different stages of the evolution 
can be significant for the evolution of the star and allow insight into the behavior of matter at 
extreme densities. 
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Neutron stars that have both hadronic and quark phases are normally described by two 
different equations of state. They are connected at the point in which the pressure of the 
quark phase becomes higher than the pressure of the hadronic phase. In our approach we 
have instead a combined model of hadrons and quarks and one equation of state for both 
phases. With this, we study important features of this transition, like the restoration of 
chiral symmetry. For this purpose we extend the generalized hadronic non-linear chiral 
model [1|2] to include quark degrees of freedom. Although the temperature in late-stage 
neutron stars is smaller than 1 MeV, in proto-neutron stars, right after the original 
supernova explosion, it can reach values up to 50 MeV. For this reason we examine the 
behavior of our model at all possible densities and temperatures. 

The lagrangian density of the non-linear sigma model, evaluated in mean field approx- 
imation, reads: 

L = LiCin + Li nt + Lsb + Lself ~ U, (1) 

where besides the kinetic energy term for baryons, quarks and leptons (included to insure 
charge neutrality), the terms: 
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Fig. 1. Left: Order parameters for chiral symmetry restoration (<r) and for deconfinement (<t>) for charge 
neutral and beta equilibrium matter at zero temperature. Right: Population for charge neutral and beta 
equilibrium matter at zero temperature 
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represent the interactions between baryons or quarks and vector and scalar mesons, 
an explicitly chiral symmetry breaking term (responsible for producing the masses of the 
pseudo-scalar mesons), the self interactions of scalar and vector mesons and a potential 
term U for the Polyakov loop as defined below. The index i denotes nucleons or up and 
down quarks. Note that although we only consider non-strange baryons and quarks in 
this study, the model includes couplings of the nucleons to the strange scalar condensate 
field C, which is taken into account here. The effective masses of baryons and quarks are 
generated by the scalar fields except for a small explicit mass term (equal to 5m = 1 50 
MeV for baryons and Sm = for quarks) and the term containing the Polyakov field $: 

ml = g ba a + g K ( + 5m b + g b ^ 2 , m* = g qa a + g q( C + Sm q + g g <s,(l - $). (3) 

As can be seen from the formula above, there are two order parameters in our model. 
The scalar condensate (a) decreases with the increase of temperature/chemical poten- 
tial reducing the effective masses and signaling the chiral symmetry restoration. The 
Polyakov loop ($) increases with increasing temperature/chemical potential, which in 
turn generates higher (lower) effective masses of the baryons (quarks), effectively de- 
scribing deconfinement. This behavior can be seen for neutron star matter in FigfT] The 
Polyakov loop is related to the Z(3) symmetry, that is spontaneously broken by the pres- 
ence of quarks. It has the form $ = jrTr[exp (i J drA^)), where A4 = iAo is the temporal 
component of the SU(3) gauge field. 

The potential for the Polyakov loop U reads: 

U = (a a T 4 + ai [i 4 + a 2 T 2 ^ 2 )^ 2 + a 3 T 4 In (1 - 6$ 2 + 8$ 3 - 3$ 4 ), (4) 
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Fig. 2. Left: Entropy profile for fixed temperatures for charge neutral and beta equilibrium matter. Right: 
Temperature profile for fixed entropies for charge neutral and beta equilibrium matter 

that was based on [31415] and adapted to include also terms that depend on the chemical 
potential. The coupling constants for the hadrons are chosen to reproduce the vacuum 
masses of the baryons and mesons and the nuclear saturation properties. The vacuum 
expectation values of the scalar mesons are constrained by reproducing the pion and kaon 
decay constants. The coupling constants for the quarks are chosen to reproduce lattice 
data and known information about the phase diagram (a first order phase transition for 
pure gauge at T — T c and a critical point dividing a crossover area from a first order 
phase transition line are successfully described in this ansatz [B]). 

An interesting effect of the deconfinement is that because of the increase of degrees 
of freedom, there is a jump on the entropy at a certain chemical potential for a fixed 
temperature (Fig. [2] left panel). As a consequence of this, if we fix entropy instead, as 
normally done in neutron star calculations [2], there will be a jump in the temperature 
(Fig.[2]right panel). This calculations were done requiring charge neutrality independently 
in each phase. If one considers global neutrality instead, the effect of the deconfinement 
transition will be smoothed out. 

Until now, hyperons were not considered because they appear at about the same 
densities as quarks and would be suppressed by the new phase. The strange quark, 
having a much higher mass, was also not considered because it only appears at very high 
densities and would practically not affect neutron star properties. More precise results 
including both these features and a full analysis about star stability will be shown in 
future work [6]. 
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